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Abstract 

We propose a computational scheme for the ab initio calculation of Wannier functions (WFs) for 
correlated electronic materials. The full-orbital Hamiltonian H is projected into the WF subspace 
defined by the physically most relevant partially filled bands. The Hamiltonian H WF obtained 
in this way, with interaction parameters calculated by constrained LDA for the Wannier orbitals, 
is used as an ab initio setup of the correlation problem, which can then be solved by many-body 
techniques, e.g., dynamical mean-field theory (DMFT). In such calculations the self-energy operator 
is defined in WF basis which then can be converted back into the full-orbital Hilbert space to 
compute the full-orbital interacting Green function G(r, r', e). Using G(r, r', e) one can evaluate the 
charge density, modified by correlations, together with a new set of WFs, thus defining a fully self- 
consistent scheme. The Green function can also be used for the calculation of spectral, magnetic 
and electronic properties of the system. Here we report the results obtained with this method 
for SrV03 and V2O3. Comparisons are made with previous results obtained by the LDA+DMFT 
approach where the LDA DOS was used as input, and with new bulk-sensitive experimental spectra. 

PACS numbers: 78.70.Dm, 71.25.Tn 
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I. INTRODUCTION 



Model Hamiltonians used in the study of correlation effects in solids have a Coulomb 
interaction term in a site-centered atomic-like orbital basis set which is not explicitly defined. 
When the correlated electrons are well localized, as, for example, 4/-states of rare-earth 
ions, atomic orbitals (or atomic sphere solutions like muffin-tin orbitals) are a good choice. 
However, the most interesting problems occur in the regime of metal-insulator transitions, 
where the states of interest become partially itinerant and rather extended. The error of 
using atomic orbitals is most severe in the case of materials with strong covalency effects, 
like late transition metal oxides, where partially filled bands are formed by the mixture of 
metallic rf-orbitals and oxygen p-orbitals. For example, in high-T c cuprates correlated states 
have the symmetry of Cu-3d x 2 — y 2 orbitals, but are actually Zhang- Rice singlets formed by 
the combination of oxygen p-states centered around the Cu ion and having x 2 —y 2 symmetry. 

In model calculations the problem of defining the correlated orbitals is not very important, 
because it only affects model parameter values, which in any case are considered fitting 
parameters. However any attempt to construct an " ab initio" calculation scheme requires 
an explicit definition of the basis set for the Coulomb interaction term. An important 
requirement for such a choice is that the orbitals must produce the partially filled bands 
where Coulomb correlations occur while preserving the localized, site-centered atomic-like 
form. These requirements are fulfilled for Wannier functions (WFs) \W^) defined as a 
Fourier transformation of the Bloch functions \ipnk) jjj. Here and below functions are labeled 
with band index n, lattice translation vector T and wave vector k. 

When there is more than one band crossing the Fermi level, WFs are not uniquely defined. 
Any k-dependent unitary transformation of the set of Bloch functions \ipnk) for these 
bands produces a new set which can be used for the calculation of WFs via Fourier trans- 

'. f one imposes the requirement that the WFs should have the 
DD , this unitary transformation is well defined. The explicit 
form of the WFs allows one to compute Coulomb interaction parameters in constrained local 
density approximation (LDA) calculations. 

In this way the parameters for the ab initio many-body Hamiltonian (non-interacting 
Hamiltonian H WF and Coulomb interaction) in the WF basis can be computed by any 
first principle electronic structure calculation scheme (below we use the LMTO method). 



formation (eq. Sec. Ill A|) . 
symmetry of atomic orbitals 
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This Hamiltonian can then be further investigated by one of the methods developed in 
the many-body community In the present work we use the dynamical mean-field theory 
(DMFT) 0, Q 0|. Within DMFT, the effective impurity problem corresponding to the 
many-body Hamiltonian is solved by quantum Monte-Carlo simulations (QMC) |7|. The 
DMFT part of the proposed calculation scheme is essentially the same as the one used in the 
recently developed LDA+DMFT approach for the ab initio investigations of correlated 
electron materials jjj. However, here we propose a more general procedure to compute the 
Green function using the Hamiltonian matrix and an integral over Brillouin zone instead of 
the Hilbert transform of the LDA density of states (DOS). This particular method allows 
one to avoid the uncontrollable errors occuring in the computation of the Green function 
using the Hilbert transform of the LDA DOS. Thus, to obtain an insulating solution we 
need to cut off the long (metal-oxygen) hybridization tails of the DOS, renormalize it and 
shift the Fermi energy to get an integer filling. In the present method we overcome the 
above-mentioned difficulties owing to the integer filling of Wannier orbitals. The result of 
the DMFT calculations is a local self-energy operator E(e). In our scheme this operator 
is defined in the WF basis set {W n } and is acting in the subspace of partially filled bands 
which are used for the construction of the WFs. 

The paper is structured as follows. In Sec. |H] the details of our scheme are presented. In 
Sec. Ill Al we describe the construction of WFs, as well as the ab-initio Hamiltonian matrix 
within this basis set in terms of Bloch functions. In Sec. IIIBI we propose a general method 
for the construction of WFs using the Green function G(r, r', e), which reduces to the results 
of Sec. Ill Al in the non-interacting case. The reason for doing so is that the correlation effects 
can significantly renormalize the electronic states of the partially filled bands. Hence the 
WFs computed from non-interacting Bloch states are not an optimal choice for the basis set 
any more. In Sec. Ill CI we discuss how to calculate within DMFT the local Green function 
with the input of the Hamiltonian matrix in WF basis set instead of the LDA DOS (which 
is valid only in the case of degenerate bands). In Sec. Ill Dl we show that the self-energy 
operator within the WF subspace, which is the solution of the correlation problem, can be 
transformed back into the full-orbital Hilbert space, thus enabling the computation of the 
full interacting Green function G(r,r',e). It can then be used to calculate the spectral, 
magnetic and electronic properties of the system under investigation. In addition, to make 
the calculation scheme fully self-consistent, one can employ the G(r,r',e) to calculate the 
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correlation affected charge density and thus the new LDA potential. Thereby the feedback 
from DMFT to LDA can be incorporated in a well-defined way. This is actually one of the 
great advantages of using the WF basis since in the LMTO basis the feedback from DMFT 
to LDA is essentially uncontrolled. In Sec. Illll the results for the electronic structure of the 
two vanadium oxides SrVC>3 and V2O3 obtained by the method developed in this work are 
presented and compared with the previous calculations by the simpler methods and new 
bulk-sensitive spectra. Finally in Sec. II VI we close this work with a conclusion. 



II. METHOD 



Let us consider the general case of the electronic structure problem. For the LDA Hamil- 
tonian H we have a Hilbert space of eigenfunctions (Bloch states \ipik)) with the basis set 
\<f>n) defined by particular methods (e.g., LMTO jlfll ]. or linearized augmented plane waves 



(LAPW) etc.). In this basis set the Hamiltonian operator is defined as 

H = J2\^) H ^(U (1) 

Here and later greek indices are used for full-orbital matrices. 

If we consider a certain subset of the Hamiltonian eigenfunctions, for example Bloch states 
of partially filled bands |^„k), we can define a corresponding subspace in the total Hilbert 
space. The Hamiltonian matrix is diagonal in the Bloch states basis. However, physically 
more appealing is a basis set which has the form of site-centered atomic orbitals. That is 
a set of WFs \W^) defined as the Fourier transformation of a certain linear combination of 
Bloch functions belonging to this subspace (see below ©). The Hamiltonian operator H 
defined in this basis set is 

H WF = J2\ W n) H nAT)(W^\. (2) 
nn'T 

The total Hilbert space can be divided into a direct sum of the above introduced subspace 
(of partially filled Bloch states) and the subspace formed by all other states orthogonal to 
it. Those two subspaces are decoupled since they are the eigenfunctions corresponding to 
different eigenvalues. The Hamiltonian matrix in the WF basis (i.e., a collection of the bases 
of the specific subspaces) is block-diagonal so that the matrix elements between different 
subspaces are zero. The block matrix H nn > in (J2J) corresponding to the partially filled bands 



can be considered as a projection of the full-orbital Hamiltonian operator JT} onto the 
subspace defined by its WFs. 

All this concerns the non-interacting (or LDA) Hamiltonian. To treat Coulomb correla- 
tions one also needs a definition of the localized orbitals where the electrons interact. WFs 
are a natural choice for such a definition. This choice leads to an important flexibility in 
the size of the basis set in the sense that the number of WFs can be changed by changing 
the set of Bloch bands considered. The simplest case is a set of partially filled bands, for 
example the ^g-bands of vanadium oxides. This is a physically justified approximation be- 
cause the Coulomb interaction happens mainly between electrons (or holes) in the partially 
filled bands. If the problem to be solved concentrates on the excitation spectrum in a small 
energy window around the Fermi level, this basis set is sufficient. However, if the excitations 
to higher lying states (real or virtual) are also important, the set of Bloch bands used to 
construct the WFs need to be extended so that the Coulomb interaction will be treated in 
a larger Hilbert subspace. 

Practically this means that the correlation problem is solved using a non-interacting few- 
orbital Hamiltonian H WF (J2J) instead of the full Hilbert space Hamiltonian H (JJJ. The 
interaction matrix elements of the model Hamiltonian can be determined from constrained 
LDA calculations for the specific WF basis set (1251) . 

Projecting the full orbital Hilbert space Hamiltonian H ((TJ onto the subspace of the par- 
tially filled bands gives us a few-orbital Hamiltonian H WF (|2|). This significantly decreases 
the complexity of the correlation problem, thus permitting its explicit solution. The many- 
body problem with a Hubbard interaction then leads to a local self-energy operator T, WF (e) 
which is naturaly defined in the basis of WFs centered on the same site: 

fWF {e) = J2\W°)Z n Ae)(W°,\. (3) 

nn' 

We note that, in contrast to other "basis-reducing" methods, the information about the 
states corresponding to the bands below and above the projected ones is not lost. In fact, 
the information is stored in the k-dependent projection matrix between the full orbital 
basis set and the orthonormalized WFs (|16J1 . The definition (J3J) allows one to convert the 
self-energy matrix H nn i[e) back to the full Hilbert space basis set (subsection III D|l . With 
this the interacting Green function can also be calculated in the full-orbital Hilbert space 
(subsection III Bj) . 
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A. Definition and construction of Wannier functions 



The concept of WFs has a very important place in the electron theory in solids since 
its first introduction in 1937 by Wannier Q. WFs are the Fourier transformation of Bloch 
states j^ik) 

\w7) = ^£ e " lkT l^>' ( 4 ) 



N k 

where N is the number of discrete k points in the first Brillouin zone (or, the number of cells 
in the crystal). These extremely convenient orthogonal functions were widely investigated 
in the seventies The strongly localized nature of the WFs together with all advantages 
of the atomic functions makes them a very useful tool where the atomic character of the 
electrons is highlighted. Thus, using the WF method, significant progress was achieved in 
the fields of narrow-band superconductors, disordered systems, solid surfaces, etc. Several 
methods for calculating WF for single and multiple bands in periodic crystals, and their 
generalization to non-periodic systems were proposed. The problem of non-unique definition 
of WFs in these methods was resolved by an iterative optimization of trial functions which 
have the same real and point group symmetry properties as WFs. Among these methods 

nn 

there are t 
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re variational Koster-Parzen principle |l5l . Ilfil ] which was generalized by Kohn [1 
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formalism proposed by Anderson |2J|, and the 
Hi- However, all these computational 



231 ] . the general pseudopotentia 
projection operator formalism by Cloizeaux 
schemes are restricted to simple band structures. 

Wannier functions are not uniquely defined because for a certain set of bands any orthog- 
onal linear combination of Bloch functions can be used in (j3J). In general it means that 
the freedom of choice of Wannier functions corresponds to freedom of choice of a unitary 
transformation matrix U^f 1 for corresponding Bloch functions Q]: 

*o = E^ k W>- (5) 

The resulting Bloch function |^k) will generally not be an eigenfunction of the Hamiltonian 
but has the meaning of a Bloch sum of Wannier functions (see below |W re k) in (JBJ)). There 
is no rigorous way to define U^f 1 . This calls for an additional restriction on the properties 
of WFs. Among others Marzari and Vanderbilt proposed the condition of maximum 
localization for WFs, resulting in a variational procedure to calculate . To get a good 



initial guess the authors of Q proposed choosing a set of localized trial orbitals \<j) n ) and 
projecting them onto the Bloch functions \ijjik). It was found that this starting guess is 
usually quite good. This fact later led to the simplified calculating scheme proposed in 
[31] where the variational procedure was abandoned and the result of the aforementioned 
projection was considered as the final step. The approach of |2J has recently been used for 
the investigation of the row of 3d transition metals (Fe, Co, Ni, and Cu) within the simplest 
many-body approximation, namely the unscreened Hartree-Fock approximation 28| . 

Another possibility to construct WFs was recenly developed by Andersen et al. (29 1. They 
proposed the Nth-order muffin-tin orbital (NMTO) scheme in which Wannier-like low-energy 
MTOs can be designed a priori. Using a new implementation of the LDA+DMFT approach 
they performed an investigation of the Mott transition in orthorhombic 3d 1 perovskites j^] ■ 
In this approach a realistic Hamiltonian constructed with Wannier orbitals (on symmetrically 
orthonormalized NMTOs) was solved by DMFT, including the non-diagonal part of the on- 
site self-energy. 

Our projection procedure works as follows. First of all one needs to identify the physically 
relevant bands which will then be projected onto a WF basis. For example, in perovskites 
one usually takes the partially filled d-shell or some particular <i-bands of transition metals, 
since they are mainly responsible for the physical properties of the system Qj. These orbitals 
are well-separated and are, in our approach, easily extracted from the full orbital space as 
will be shown later. Moreover, the projection method is applicable even in the case where 
the bands of interest differ and are strongly hybridized (for example, Cu-3d and 0-2p states 
in high-T c superconductors jsij). 

To project bands of particular symmetry onto the WFs basis one can select either the 
band indices of the corresponding Bloch functions (Ai, A2), or choose the energy interval 
(Ei, E 2 ) in which the bands are located. Non-orthogonalized WFs in reciprocal space |W n k) 
are then the projection of the set of site-centered atomic- like trial orbitals |0 n ) on the Bloch 
functions \ipik) of the chosen bands (band indices Ai to A 2 , energy interval (_Ei,_E 2 )): 

N 2 

iWnk) = ^2 \1>ik)W*\<t>n) = Yl l^k)(^k|0n). (6) 

i=Ni »(£a<ei(k)<^2) 

Then the real space WFs \W^) are given by 

i^ T ) = ^E e ~ ikT i^>- ( 7 ) 



N k 



s 



In the present work the trial orbitals \<p n ) are LMTOs. Note that in the multi-band case 
a WF in reciprocal space |W n k) does not coincide with the Bloch function |^ n k) due to 
the summation over band index i in (jHJ). One can consider them as Bloch sums of WFs 
analogous to the basis function Bloch sums <p^(r) ©• 

The coefficients (ipik\4> n ) in (jHJ) define (after orthonormalization) the unitary transforma- 
tion matrix in @. However, the projection procedure defined in (JHJ) is more general 
than the unitary transformation (0). Namely, the number of bands (N2 — N\ + 1) can be 
larger than the number of trial functions. In this case the projection © will produce N new 
functions |W„k) which define a certain subspace of the original (N2 — N± + l)-dimensional 
space. This subspace will have the symmetry of the set of trial functions. In the next 
subsection we propose a way to determine WFs from the Green function of the system (J28|) 
rather than from a set of Bloch states as in ©. In this alternative projection procedure, 
trial functions are projected onto the subspace defined by the Green function in a certain 
energy interval. 

The Bloch functions in LMTO basis (or any other atomic orbital-like basis set) are defined 

as 

\m = J2° k M)> (8) 

where \i is the combined index representing qlm (q is the atomic number in the unit cell, Im 
are orbital and magnetic quantum numbers), 0„(r) are the Bloch sums of the basis orbitals 

M* ~ T ) 

* r ) = ^E^^-T), (9) 
and the coefficients have the property 

c% = <0 M |Vik). (10) 

If n in \<p n ) corresponds to the particular qlm combination (in other words \<p n ) is an 
orthogonal LMTO basis set orbital), then (^ik|0n) = c^*, and hence 

N 2 N 2 

\w nk ) = ^ i^ ik >4* = x; E4*l^> = £%»W>> ( n ) 

i=Ni i=Ni n /j, 

with 

%n = E ^nt- (12) 

i=Ni 



For a non-orthogonal basis set see Appendix O 

In order to orthonormalize the WFs (fTTj) one needs to calculate the overlap matrix O nn >(k) 

N 2 

O nn ,(k) = (W nk |^ k ) = C ni^, (13) 

i=Ni 

and its inverse square root S nn r(k) is defined as 

Snn'(k) = 0_Y 2 (k). (14) 

In the derivation of (fT3*j) the orthogonality of Bloch states (^nklVVk) = Snn' was used. 
From (fTTj) and the orthonormalized WFs in k-space |W rek ) can be obtained as 

\W nk ) = J2SnnW\W n , k ) = 2 |^k>4* = Y. h lnK)> (15) 



i=JVi 



with 



-k* 



• =(^k|^nk) = ^^(k) C ^, (16) 



N 2 

= (^\W nk ) = £ <*c£. (17) 

i=JVi 

The real space site-centered WFs at the origin \W®) are given by the Fourier transform 
of |W nk ) with T = 0. From (JljJ) and © one finds 

ViV k T> V k 7 

= 5>'(n,/i,T)0 M (r-T) (18) 
= ^w(n,s)0 a(s) (r - T s ), 



where w' and u> are the expansion coefficients of WF in terms of the corresponding LMTO 
orbitals, in particular, 

^M^E^^K (I 9 ) 

k 

Here s is an index counting the orbitals of the neighboring cluster for the atom where orbital 
n is centered (T s is the corresponding translation vector, a(s) is a combined qlm index). 
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The explicit form of the real space WF ()18|) can be used to produce, e.g., shapes of chemical 
bonds. 

For other applications only the matrix elements of the various operators in the basis of 
WF ()15jl are needed. From (|15|) . and (J 18)) the matrix elements of the Hamiltonian H WF 
in the basis of WF in real space where both orbitals are in the same unit cell are 

/ 1 N2 \ 
HZf(0) = (^lbEEl^(k)W 

^ k i=iVi ' 
1 ^2 

= jyEE^( k )- ( 2 °) 

k i=Ni 

6j(k) is the eigenvalue for a particular band. 

If, on the other hand, one of the orbitals corresponds to the WF for the atom n' shifted 
from its position in the primary unit cell by a translation vector T, then the corresponding 
Hamiltonian matrix element is 

, ^2 

HZ?(T) = (W2\H\W%) = - E E *^ ei (k)e- ikT . (21) 

k i=Ni 

Matrix elements of the density matrix operator (occupation matrix Qj^f ) in the basis of 
WFs can be calculated as 

Qinv = (wz\ T7 E E - c(k)x^i \w%) = 



N 

k i=Ni 



„ E E **W - Q(k))e- kT (22) 



k i=JVi 



6{x) is step function, Ef is the Fermi energy. 

Finally, the matrix elements of the Hamiltonian H WF in reciprocal space are 

^(k) = <W„ k | ^EEl^)^(k')(^k'| |W„, k ) 

^ k' i=iVi ' 

= E^e^k). (23) 

t=JVi 

The is valid only if the WFs computed by eqs. (|T3j) - (fT7|) . If the WFs were obtained 
in one calculation and then used to compute the Hamiltonian matrix in another (as is the 
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case for the WFs ()28|) in the Green functions formalism (see subsection IIIB|) ) then eq. 
is not valid any more and the general expression must be used: 

HifiX) = E e *( k ) E b >% E ( 24 ) 

i=Ni (J, v 

Thus, the transformation from LMTO to WF basis set is defined by the explicit form of 
WFs (J15I17J) . and by the expressions for matrix elements of the Hamiltonian and density 
matrix operators in WF basis ()2(JI22)) . The back transformation from WF to LMTO basis 
can also be defined using ()15|) (see subsection HID)) . 

Finally, the Coulomb matrix element U needs to be calculated in the same WF basis. 



This requires a method similar to constrained LDA [31|, but now for WFs. To this end the 
WF-energy (|20|) is computed as a function of its occupancy ()22|) for a given WF n. Then 
the corresponding Coulomb interaction parameter U n in the WF basis is given by 

rr = dHZ F (0) 

n ~dQZ F (oy (25) 

As one can see U n depends on the WFs via (|20I22|) . Once the WFs have been recalculated 
(for example in some self-consistent loop) the interaction has to be recalculated as well. 

B. Wannier functions in the Green function formalism 

In many-body theory the system is usually not described by Bloch functions © 
and their energies ej(k) but by the Green function 

G(v, r', e) = G k (r, r', s) = ^ E €^K^)^')- ( 26 ) 

k k iiv 

The Green function matrix Grt v [e] is defined via the non-interacting Hamiltonian matrix 
H^ik) and the self-energy matrix Ej^(e) as 

G^M = (e- H(k) - E(e, k) + zr^ 1 . (27) 

We define non-orthonormalized WF obtained by projecting the trial orbital n (r) on the 
Hilbert subspace defined by the Green function (|26[) in the energy interval (Ei, E 2 ), namely, 

^nk(r) = -hm J de j dr'G k (r, r', e)^(r') = E?n^W, (28) 
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and 



%. = ~Im I deG* n (e). (29) 



Ex 



In the non-interacting case, the self-energy operator E(e, k) is absent, and hence we have 

^-Et^Ht^ < 3 °> 

where 4 are the eigenvectors (fTOj). and e$(k) are the eigenvalues of -f^(k). Thus b^ n in (f2"9l 

becomes 

^2 

bl n = £ 44*, (31) 

i=JVi 

where iVi, A2 are the band numbers which correspond to the energy interval (E±, E2). Since 
this recovers the result of (fT2|) . we demonstrated that our general definition of WFs (J2Hj) via 
Green function reduces to that in terms of Bloch functions ()11|) in Sec. Ill Al 

To orthonormalize W n k(r) defined in ()28|) . one can just follow the orthonormalizing pro- 
cedure made in Sec. Ill Al (I13H17|) . which will not be repeated here. But it should be pointed 
out that in the Green function formalism the overlap matrix O nn i (k) is defined as 

The occupancy matrix in the orthogonalized WF basis (f2*K|) is defined as 

Ep 

Qn,AT) = -Urn J dej j rfrdr'l^W r : k (r)G' k (r,r / ,£)iy n , k (r')e- ikT . (32) 
-00 k 

By using (|2T)j) and orthogonalized (|2~%)). one finds 

Qnn' (T) b >vn>Q>~ lilT , (33) 



with 



N 

k \iv 



Ql„ = -Um J deGlM (34) 

— oo 
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The energy matrix can be defined similarly (except that the integral over energy is cal- 
culated in the {E^E-i) interval where the corresponding WFs are defined) as: 



E nn r ( T ) 



l -Im J ede J J dvdr'^ ^ WC k (r)G k (r, r', e)WV k (r')e-* kT (35) 

Ei k 



-ikT 

k fiu 



with 



E 2 

1 



E%, = —Im / eefeGjUe)- (36) 



7T 

Si 

While ()35|) looks similar to the non- interacting Hamiltonian in WFs basis (J2T|) . it includes 
correlations via S(e) in (|27|) and hence is interacting. 

C. DMFT in the Wannier function formalism 

In the previous subsection we showed (eqs. EoTI^Hj) that the self-energy operator is needed 
to construct the WFs in terms of the full interacting Green function. The DMFT 0,0,1(1 was 
recently found to be a powerful tool to numerically solve multi-band Hubbard models. To de- 
fine parameters of the correlated model Hamiltonian (hoppings, screened Coulomb integrals), 
density functional theory within the LDA was used sj]. The combined LDA+DMFT com- 
putational scheme was successfully applied to a wide range of compounds with degenerate 
(or almost degenerate) orbitals (for more details see Q]). In these cases the non-interacting 
LDA DOS was used to obtain the Green function of the system through a Hilbert transfor- 
mation. Furthermore, the screened Coulomb interaction parameters U and J were calculated 
by constrained LDA j3jj. 

Quite generaly, this scheme needs to be improved in two respects: (i) instead of the LDA 
DOS an LDA-Hamiltonian with a few, relevant orbitals should be used to calculate the 
Green function, and (ii) a feedback from DMFT to LDA should be incorporated. Both of 
these problems are solved by the new approach proposed in this work. In this method the 
Hamiltonian matrix in the WF basis set n^J (k) is calculated from the LDA Hamiltonian 
via the projection procedure (jHESJ)- In the DMFT self-consistency loop the local Green 



14 



function G °^,(e) is then computed as an integral over the first Brillouin zone (BZ) : 



G 



loc 



1 



V, 



BZ 



dk 



H WF (k) - E WF (e 



WF, 



-1 



(37) 



The integration can actually be restricted to the irreducible part of the BZ via the analytical 
tetrahedron method j^j with a subsequent symmetrization of the Green function matrix. 
The chemical potential E^ is determined by the number of electrons on the N interacting 
orbitals of interest 3^ |. 

The DMFT is based on the fact that in the d = oo limit the self-energy operator is 
local (34, |2| . Its matrix S n „/(e) (n, n' - WF indices) is defined in WF basis (J3J). If the trial 
functions in (|6I28|) are chosen as the basis functions of the irreducible representation of the 
point symmetry group of some particular real system j^, the Green function matrix (|27J1 
and hence the self-energy matrix (J3J) can be made diagonal [3^ in the n index for on-site 
matrix elements. 

To set up the DMFT equations one needs to define the bath Green function Q in the 
usual way via the Dyson equation 

g-l = (Qlocyl + S> ( 3g) 



To obtain the Green function G imp of the effective single impurity Anderson problem, vari- 
ous methods can be used: quantum Monte-Carlo (QMC), numerical renormalization group 
(NRG), exact diagonalization (ED), non-crossing approximation (NCA), etc. (for a brief 
overview of the methods see j|). With the condition G loc = G tmp one closes the self- 
consistent loop which can then be iterated until a converged solution for the self-energy 
Ennfe) is found. 



D. Converting back to the full-orbital Hilbert space 

The self-energy operator T, WF (e) obtained as a solution of DMFT in Sec III CI is defined in 
the WF basis set (J5J). In order to compute the interacting Green function in the full-orbital 
Hilbert space (|2fill27j) one has to convert it back to the full-orbital (LMTO) basis set. This 
can be easily done by using the linear expansion form of the WFs in terms of the full-orbital 
basis set (jl5H7j) . 

S£,(e) = (#|S(e)|#) = ^(0^k>£nn^)<^»'k|0*> = ]T^ n £ nr ^)C- (39) 

n n 
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Here we use the local form of the self-energy operator as obtained in DMFT, but the for- 
malism can be easily generalized. In the following we refer to this self-energy operator as 
the "full-orbital" self-energy. 

The matrix elements of the self-energy operator T^ u (e) (|39|) together with the non- 
interacting Hamiltonian matrix allow one to calculate the Green function matrix G* u (e) 
(|2*7)l and thus the full-orbital interacting Green function G(t,t',s) (J2SJ)- G(r,r',e) contains 
the full information about the system, and various electronic, magnetic and spectral proper- 
ties can be obtained from it. In Sec. IIIIBI we use the full-orbital interacting Green function 
computed within DMFT(QMC) to calculate the photoemission and X-ray absorption spec- 
tra for the strongly correlated vanadium oxides SrVOs and V2O3, and to compare them with 
new bulk-sensitive experimental spectra. 

One can also calculate the charge density distribution modified by correlation effects via 

p(r) = ~Im J deG{r,r,e). (40) 

—00 

With this p(r) one can recalculate the LDA-potential (which is a functional of electron 
density). From the full-orbital Green function (ffilj) one can recalculate new WFs ()28I29|) 
which together with the new LDA Hamiltonian allows one to obtain new parameters for 
the non-interacting Hamiltonian (|24|). With (|2*3j) one can then compute a new Coulomb 
interaction parameter U. The set of new LDA potential, WFs and Coulomb interaction 
parameters calculated from the interacting Green function (J26|) defines the input for the 
next iteration step and hence closes the self-consistency loop in the proposed computation 
scheme. For the feedback from DMFT to LDA in the particular case of the LMTO method 
[lol | one needs a set of moments for the partial densities of states for every atomic 

sphere q and the orbital moment / in order to calculate the new charge density and 
hence the new LDA potential: 

Ep 

M< m) = / dse m N ql (e), (41) 
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FIG. 1: Scheme of a& initio fully self-consistent LDA+DMFT scheme based on the WF formalism 
(see text). Red color marks steps that were not performed in this paper. 

E. Summary of the WF scheme 

For clarity, in Fig. ^the essential steps of the WF scheme presented here are summarized. 
There are four interconnected parts in this scheme: (i) the basis of WF, (ii) the matrix 
elements of the Hamiltonian and the self-energy in WF basis, (iii) the Coulomb interaction 
between electrons on the WFs, and (iv) the projection into the few-orbital basis and back 
transformation to the full-orbital basis which retains the information about all orbitals. 
First the matrix elements of the non-interacting Hamiltonian in reciprocal space H^f (k) 
(|23jl and the interaction term Y, nn /(iu)) (JSJ are written in the basis of explicitly defined 
WFs |Wnk) (|15jl . The actual correlation problem, defined by the sum of these two terms 
P7jl. is then solved within the LDA+DMFT(QMC) approach |gj. The local self-energy 
H nn i(ioS) obtained thereby is then transformed back from the Wannier basis to the full- 
orbital space (see subsection III D|) . Furthermore, with the full-orbital self-energy Yi^iiu) 
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(J!!!)) the full-orbital Green function G^ v {iuj) (126127)) for the correlated electrons is calculated 
by a k-integration over the Brillouin zone 3|| . The Green function allows one to determine 
a new charge density distribution (|4U|) (in the LMTO case see (|41|) ) and a new set of WFs 
(via (|28l29j) ) . This is used to construct a new LDA-potential and new non- interacting 
Hamiltonian. Together with the new Coulomb interaction parameters U n ()25|) they serve as 
the input for the next iteration of calculations, thus completing the self-consistency loop. It 
should be stressed that in this scheme self-consistency involves not only the self-energy but 
also the basis of WFs in which it is defined, the charge density and LDA potential used for 
constructing the non-interacting Hamiltonian, and the interaction strength between electrons 
in the WFs. 

After convergency is reached the maximum entropy method (MEM) j^] can be used to 
obtain spectral functions. Then, using the Kramers-Kronig transform, the Green function 
on real axis G nn i{uj) is computed. With that one can construct the WFs basis self-energy 
Ti nn i{u) (see Appendix IB2|) . By the back transformation of £ nri '(u;) with a subsequent 
k-integration the full-orbital Green function G^ v {uj) is obtained, which now contains infor- 
mation not only about the states for which correlations are considered, but also for all other 
orbitals of the system. G^(lj) is used to obtain orbitally resolved spectral functions. This 
allows one, for example, to investigate the influence of the correlated orbitals (e.g., partially 
filled t2 9 orbitals of V in SrVOs and V2O3) on other orbitals (oxygen 2p and Vanadium e g 
orbitals). It also makes possible the computation of spectral functions in a wide energy 
region and not only in the vicinity of the Fermi level. 



III. RESULTS AND DISCUSSION 



In our earlier LDA+DMFT calculation scheme the specific properties of a material en- 
tered only via the LDA partial densities of states (DOSs) for the orbitals of interest. This 
procedure is valid for systems where the bands of interest are degenerate (as in cubic crystals; 
see Appendix El . For more complicated systems with lower symmetry one needs to employ 
the scheme proposed in this work, where the non- interacting Hamiltonian H WF (k) (|2*H|) 
(projected on WFs) describing the N orbitals under consideration is used for the calculation 
of G l ™,( £ ) (jHZD within DMFT (see section ITTT1. 

In this section we present results of LDA+DMFT calculations using the projected Hamil- 
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tonian H WF (k). The scheme was applied to SrVC>3 which has a cubic perovskite crystal 
structure, and to the trigonally distorted V2O3 (both in the insulating and metallic phase). 
The results are compared with previous LDA+DMFT calculations where the Hilbert trans- 
formation of the LDA DOSs was used for SrVC>3 4(]] (see section IHI Aj) . One should note 
that the DOSs used for DMFT calculations of V2O3 were obtained by the TB-LMTO-ASA 
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code v. 47 in contrast to where the DOSs were calculated by the ASW method 
However, the DOSs obtained in both methods are very similar and do not produce much 
different LDA+DMFT results. 

The full-orbital calculation scheme proposed in this work allows one to answer an impor- 
tant question: how do Coulomb correlations between some orbitals affect the other orbitals, 
and in particular, how does the interaction between the partially filled t2 g orbitals of the 
V3<i-shell influence the filled 0-2p and the unoccupied V-3d(e°) bands? To answer these 
questions the matrix of the self-energy operator S nn /(e) Q was converted back to the full- 
orbital basis set from the WF basis, and the full orbital interacting Green function (fTTj) 
was calculated. ^From that, total and partial DOSs were computed to produce theoretical 
spectra for comparison with the experimental photoemission and X-ray absorption data. 

In this work, QMC simulations []| were used to solve the effective single impurity An- 
derson problem in the DMFT loop. The result of the DMFT(QMC) calculation is the 
self-energy on imaginary energy axis S nn /(zcj). To find the full-orbital self-energy T^ u (e) 
on the real energy axis one has to perform an analytical continuation. The procedure is 
described in Appendix El A self-consistent computation of the charge density as described 
in section HlDl was not yet performed. Investigations of correlation effects on the unoccupied 
cubic e a g states for both SrV03 and V2O3 are also a matter of future calculations. 



A. Comparison of DMFT results obtained by Hilbert transformation of the LDA 
DOS and the projected Hamiltonian H WF (k) 

The cubic perovskite SrV03 can serve as a test case for our new method. For the cubic 
Oh point symmetry group, the three 3<i-orbitals xy,xz,yz transform according to the triply 
degenerate t2 9 irreducible representation of this group. Hence the corresponding Green 
function and self-energy matrices have diagonal form with equal diagonal elements. As 
shown in Appendix [XJ the results of the Hilbert transformation using LDA DOS must 
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FIG. 2: Comparison of V-3d(t 2ff ) spectral functions for SrV0 3 calculated via LDA+DMFT(QMC) 
using: LDA DOS - light line, projected Hamiltonian - black line. Fermi level corresponds to zero. 

coincide with the results of the procedure of integration over the Brillouin zone with the 
projected 3x3 t 2fl Hamiltonian. 

We use the same LDA DOS and interaction parameters as in our previous papers 
The V-3d(t 2 c,) states form a partially filled band in SrV0 3 . All t 2g orbitals (xy, xz, zy) 
are equivalent due to the cubic symmetry of the lattice, so only the results for one of the 
t 2g orbitals are presented in Fig. El In this figure V-3d(t 2g ) spectral functions, calculated 
using the LDA DOS and the projected Hamiltonian are shown. It is easy to see that both 
results are almost identical. The small differences between these two curves are due to the 
MEM (43! used for the calculation of the spectral function on the real energy axis from the 
DMFT Green function. 

We also applied our Hamiltonain scheme to a more complicated system with lower sym- 
metry where one can expect deviations from the results obtained using the LDA DOS. We 
performed LDA+DMFT(QMC) calculations for the insulating and metallic phases of V 2 03 
le projected 3x3 t 2g Hamiltonian and several U values (the values of U are the same 



with t 
as in 



4l|). In contrast to SrVOa, V 2 03 has a non-cubic trigonal symmetry (space group 
R3c). Therefore we use the basis functions of the trigonal point group instead of the 
cubic Oh group. In this basis set the V-3<i-shall is split into two groups of bands. The three 
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t2g bands are located around the Fermi energy, the two degenerate e^ bands are at higher 
energies. Interesting to us are the partially filled t2 9 bands which are formed by one a\ g and 
two degenerate orbitals. 

In Fig. El we present V-3rf(t 2s ) spectral functions resulting from DOS (LMTO and ASW) 
and Hamiltonian calculations at U= 4.5eV, ave rag ed over the three t2 9 bands. The spectrum 



computed with ASW input is taken from Ref. 
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FIG. 3: Comparison of V-3d(t2 S ) spectral functions for V9O3 in the metallic phase calculated via 
LDA+DMFT(QMC) using: LDA DOS (LMTO) - light line |4l|, LDA DOS (ASW) - dashed line, 
projected Hamiltonian - black line. Fermi level corresponds to zero. 



As mentioned before, the differences between the curves calculated with ASW and LMTO 
DOS as input are small; there are only minor deviations in the peak height. In the compar- 
ison between the DOS calculations and the Hamiltonian calculation, it is interesting to note 
that for these averaged spectra, the differences are relatively small. We find the typical four 
peak structure (with lower Hubbard, quasiparticle peak and double-peaked upper Hubbard 
band split by Hund's rule coupling) for all three calculations. Only the position and the 
height of the quasiparticle peak are quite different for the Hamiltonian calculation, indi- 
cating a more insulating solution than for the DOS calculations at the same U-value. The 
similarity between the results for the DOS and Hamiltonian calculations is not surprising, 
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because the trigonal distortion in V2O3 is relatively small. Therefore the center of gravity 
and the bandwidth of the three t2 9 bands do not change much and the DOS calculations 
can still produce accurate results. 

The differences between DOS and Hamilton calculations are more pronounced when one 
compares the band-resolved t2 S spectra. In Fig. 0] and Fig. El the &i g and e^ spectra for the 
insulating and metallic crystal structure are presented for U= 4.5, 5.0 and 5.5 eV. 
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FIG. 4: Comparison of V-3d(t2 S ) spectral functions for V2O3 in the insulating phase calculated 
via LDA+DMFT(QMC) using: LDA DOS - light line, projected Hamiltonian - black line. Upper 
figures - a± g , lower figures - orbitals. Fermi level corresponds to zero. 



Here the curves calculated by different schemes are distinctively different. This is es- 
pecially clear from the upper part of the figures EJ El corresponding to the ai g -orbital of 
the V-3d(t2 9 ) subband. This is due to the fact that the hybridization effects are better ac- 
counted for by the projected Hamiltonian, which includes not only intra-band hoppings but 
also inter-band on-site and inter-site hoppings. The latter effectively increases the bandwidth 
of the ai s -orbital. The corresponding spectral functions have a different peak structure and 
different intensities. 

Another conclusion from this series of figures is that the results obtained with the Hamil- 
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Energy (eV) 

FIG. 5: Comparison of V-3d(t2 9 ) spectral functions for V2O3 in the metal phase calculated via 
LDA+DMFT(QMC) using: LDA DOS - light line, projected Hamiltonian - black line. Upper 
figures - dig, lower figures - orbitals. Fermi level corresponds to zero. 

tonian procedure show a more insulating behavior for the same values of U than the results 
with the scheme using the LDA DOS. 

B. Results of full-orbital calculations: comparison of calculated spectra with ex- 
perimental PES and XAS data 

The full-orbital calculations scheme proposed in this work produces an interacting Green 
function G(r, r', e) (|2f)j). The knowledge of the full Green function allows us to calculate the 
spectral function not only for the states with correlations (VSdfag) orbitals in the case of 
SrV03 and V2O3), but also their effect on the lower lying occupied oxygen 2p states and 
the higher lying unoccupied V-3d(eZ) states, thus facilitating a comparison of the calculated 
and experimental spectra in a wide energy region. 

Valence band photoemission spectroscopy (PES) using high photon energies and O 
Is — > 2p x-ray absorption spectroscopy (XAS) were performed on the beamline BL25SU 
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FIG. 6: Self-energy on real energy axis for SrVC>3. Real part - full black line, imaginary part - full 
light line. Fermi level corresponds to zero. 

at the SPring-8 synchrotron radiation facility. The PES spectra were taken using a 
GAMMADATA-SCIENTA SES-200 electron energy analyzer and XAS spectra were mea- 
sured by the total electron yield. The overall energy resolution was set to 0.2 eV. The 
pressure in the analyzer chamber was about 4 x 10 -8 Pa. Single crystals of SrVC>3 were 
cooled to 20 K and clean surfaces were obtained by fracturing in situ for the PES spectra, 
and by scraping in situ for the much more bulk sensitive XAS spectra. Well-oriented single- 
crystalline V2O3 samples were cleaved in situ at a temperature near the metal-insulator 
transition, yielding clean specular surfaces. The surface cleanliness was confirmed before 
and after the spectral run. 

We start with the results for the SrVC>3 system. In Fig. the self-energy on the real 
energy axis for SrVC>3 calculated via (jB2|) is shown. The m*/m ratio is calculated via 
m*/m = 1/(1 — (<99?X(c<j) / 'did)) . This self-energy was used for the calculation of total and 
partial DOSs in the full-orbital Hilbert space (see Fig. l7l8j) . In Fig.[7[ the total spectral func- 
tions of SrV03 are presented. Differences between LDA and full-orbital LDA+DMFT(QMC) 
spectral functions mainly occur near the Fermi level. In particular, the LDA spectral func- 
tion has a more pronounced quasiparticle peak. The DOS calculated using the full-orbital 
self-energy has a three-peak structure: lower Hubbard band (suppressed by oxygen states), 
quasiparticle peak and upper Hubbard band located at about 3 eV. The origin of the up- 
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FIG. 7: Comparison of total spectral functions of SrV03 calculated via: LDA - full light line, using 
the full-orbital self-energy from LDA+DMFT(QMC) - dashed black line. Fermi level corresponds 
to zero. 




2 
Energy (eV) 



FIG. 8: Partial V-3d spectral functions for SrV03 calculated using the full-orbital self-energy from 
LDA+DMFT(QMC). Total d - full light line, t 2g - full black line, e g - dashed black line. Fermi 
level corresponds to zero. 
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FIG. 9: Comparison of photoemission spectra of SrVO"3 with spectral functions calculated via 
LDA+DMFT(QMC): taking into account only t2g - light line, using full-orbital self-energy - black 
line. Theoretical spectra are multplied with the Fermi function corresponding to 20K and broad- 
ened with a 0.2 eV Gaussian to account for the instrumental resolution. Intensities are normalized 
for quasiparticle peaks. Fermi level corresponds to zero. 

per Hubbard band becomes clear from Fig. |H| One can see that this broad peak is the 
sum of V-3d(eg) and V-3cZ(t2 S ) states. It should be noted that correlations on the V-3d(e°) 
orbitals were not explicitly included here. We find, however, that the V-3d(e°) states are 
slightly modified by the full-orbital self-energy due to the of hybridization with the cor- 
related V-3d(t2g) orbitals. The question how correlations affect the position and width of 
V-3d(eg) states directly can only be answered by employing the full 3<i-shell 5x5 projected 
Hamiltonian. 

Introducing correlations between t2 S states changes significantly the total and partial 
LDA DOSs of SrVC>3. The main modification is a transfer of spectral weight from the 
energy region near Fermi level to the lower and upper Hubbard bands, and the reduction 
of the weight of the quasiparticle peak. Our calculations find a strongly correlated but still 
metallic ground state for SrVOs. 

To compare our results with the experimental PES we calculated a weighted sum of V-3d 
and 0-2p spectral functions according to the photoemission cross section ratio 3:1, corre- 
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FIG. 10: Self-energy on real energy axis for V2O3. Real part of the self-energy from 
LDA+DMFT(QMC) - full black line, imaginary part - full light line. Fermi level corresponds 
to zero. 

sponding to an experimental photon energy 900 eV. The theoretical spectra were multiplied 
with the Fermi function corresponding to 20 K and broadened with a 0.2 eV Gaussian to 
take into account the instrumental resolution. In Fig. |§] one can see that the full-orbital 
spectra obtained in this way describe not only the quasiparticle peak, but also the peak at 
-6 eV and the shoulder at -3.5 eV of the PES spectra. Since previous LDA+DMFT(QMC) 
results were taking into account only the t 2g states they were not able to describe PES 
spectra below -2 eV. It is interesting to note that previous experimental studies did not find 
any states at -3.5 eV; only the new spectra reported here show this feature. 

The influence of correlation effects on the electronic structure will be more pronounced 
for systems close to the Mott insulator transition. For this purpose we compare the total 



and partial DOSes for V2O3 in the insulating and metallic phases calculated via LDA and 
full-orbital LDA+DMFT 46[. Since the V-3d(t2 9 ) orbitals in trigonal V2O3 split into non- 
equivalent &i g and e^ states, the self-energy will be different for these orbitals (see fig fTUjl . 

In Fig. ITT1 one can see that the introduction of correlation effects changes the total DOS 
drastically. Whereas the LDA DOS is metallic for both metallic and insulating phases, 
the LDA+DMFT spectra clearly show the metal-insulator transition. Another important 
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FIG. 11: Comparison of total spectral functions of V2O3 calculated via: LDA - full light line, using 
full-orbital self-energy from LDA+DMFT(QMC)- dashed black line. Fermi level corresponds to 
zero. 

feature are the modifications near the Fermi energy and in the energy window from -3 eV to 
7 eV. Moreover, there are also changes in the e a g states which are caused by the hybridization 
with the correlated t 2g bands, and are not due to a direct calculation of correlations in the 
e° states. In Fig. El the total d-state DOS above the Fermi energy, consisting of t 2g and e a g 
states is shown. 

In Fig. ED a comparison of experimental PES and calculated LDA+DMFT(QMC) (t 2g 
only and full-orbital) spectra is presented. The full-orbital spectrum is a weighted sum 
of the calculated V-3d and 0-2p spectral functions, according to the photoemission cross 
section ratio 2:1, corresponding to the experimental photon energy of 500 eV. The theoretical 
spectra were multiplied with the Fermi function and broadened with a Gaussian of 0.2 eV. j^J 
Below the Fermi energy the LDA+DMFT(QMC) spectral functions (for t 2g only and the 
full-orbital scheme) agree quite well with the PES spectra. However, the theory curves 
do not yet describe fine details. Specifically, the PES spectrum shows two definite slope 
changes, at roughly -0.6 eV and -1.5 eV, producing a rather flat topped spectrum centered 
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FIG. 12: Partial V-3d spectral functions for V2O3 calculated using full-orbital self-energy from 
LDA+DMFT(QMC). Total d - full light line, t 2g - full black line, e g - dashed black line. Fermi 
level corresponds to zero. 

on -1 eV, whereas the theory curves show a single peak centered on -1 eV. 

A comparison of the calculated LDA+DMFT(QMC) full-orbital spectral functions and 
the O Is — > 2p XAS spectrum is shown on Fig. The full-orbital spectra (which are the 
partial 0-2p spectra for the unoccupied V-3d states) are found to be in rather good agreement 
with the structures of the experimental spectrum their relative intensities and their energies. 
This is due to the inclusion of 0-2p and V-3d(e^) states in the calculations. The strong 
hybridization of the 0-2p with the d-states is described more correctly in the full-orbital 
calculations and the inclusion of e a g states (even without correlations) significantly changes 
the V-3<i-shell in the energy regions above the Fermi edge. The inclusion of correlations in 
e a g states within a 5x5 projected V-3d Hamiltonian is expected to add spectral weight in the 
energy interval between 3 eV and 5 eV (see Fig. ITlj) . 
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FIG. 13: Comparison of photoemission spectra of V2O3 with spectral functions calculated via 
LDA+DMFT(QMC): taking into account only t2 9 - light line, using full-orbital self-energy - black 
line. Theoretical spectra are multiplied with the Fermi function and broadened with a 0.2 eV 
Gaussian to simulate instrumental resolution. Intensities are normalized for peaks situated around 
-1 eV. Fermi level corresponds to zero. 

IV. CONCLUSION 

We formulated a fully ah initio and self-consistent computational scheme based on Wan- 
nier functions (WF) for the calculation of the electronic structure of strongly correlated 
compounds. The WF formalism provides an explicit strategy for the construction of the 
matrix elements of the required operators (Hamiltonian, self-energy operator, etc.), both in 
full-orbital and few-orbital bases, in real and reciprocal representations. The WF formalism 
allows one to project these operators from the full-orbital space to the few-orbital space and 
back, keeping the complete information about all orbitals. These projections are the essential 
ingredients of the computational scheme presented here. The self-consistency involves not 
only the self-energy but also the WF basis itself, the charge density with the LDA-potential 
and the interaction strength parameters between the electrons on the WFs. The spectra 
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FIG. 14: Comparison of O Is — ► 2p x-ray absorption spectroscopy (XAS) spectra of V2O3 with 
spectral functions calculated via LDA+DMFT(QMC) using full-orbital self-energy. Theoretical 
spectra are multiplied with the Fermi function and broadened with a 0.2 eV Gaussian to simulate 
instrumental resolution. Intensities of the theory curves are normalized to the correct number of 
V-3d electrons in the unit cell (8 electrons). The experimental curve is normalized to the same 
peak height as the full orbital curve. Fermi level corresponds to zero. 

obtained thereby found to be in good agreement with new bulk-sensitive experimental data. 

In the present work we did not yet employ the full scheme (see Fig. [IJ to investigate 
spectral functions of strongly correlated systems but employed only few-orbital Hamiltonians 
with t2 9 symmetry. Clearly, full ci-shell DMFT(QMC) results including e g states will provide 
additional information about correlation effects in the system. Such studies, as well as 
constrained LDA calculations in the WF basis and investigations of the feedback from DMFT 
to LDA part are in progress now. 
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APPENDIX A: HILBERT TRANSFORMATION 



For cubic systems the matrix of the self-energy operator (for example for t 2 g orbitals) is 
diagonal and all diagonal elements are equal. Therefore the calculation of the Green function 
within DMFT by integration of the Hamiltonian over the BZ is equivalent to the Hilbert 
transformation of the non- interacting (LDA) DOS N°(e): 
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APPENDIX B: SELF-ENERGY ON THE REAL ENERGY AXIS 



DMFT produces Green functions and self-energies on the imaginary energy axis. With 
the maximum entropy method the spectral function on the real energy axis is calculated, 
which yields the imaginary part of the Green function. In order to obtain the self-energy on 
the real energy axis, the full complex Green function G(e) is calculated using its imaginary 
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part obtained by MEM: 
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The self-energy for real energies is then calculated by solving the following two equations 
with the two variables _?__(_:) and 3E(e) : 

SR, 3{G(e)} = ^{J {e- H(k) - E^))" 1 ^}, (B2) 

BZ 

where 

5_(e) = 8___(e) + iM(e). (B3) 



APPENDIX C: NON-ORTHOGONAL BASIS SET 

Eq. is valid only for orthogonal LMTO orbitals. In the case of general non-orthogonal 
LMTOs (or any other atomic type orbital basis set), an orthogonalization procedure can be 
used by defining an orthogonal Hamiltonian H and the corresponding eigenvectors C for the 
non-orthogonal Hamiltonian H and overlapping matrix O: 

H = 0~*HO~*, (CI) 
C = 0*C. 

This orthogonalization is equivalent to the basis set transformation: 

I0n)=^oj|0 n ,)- (C2) 

n> 

Then for the non-orthogonal LMTO the trial function |0 n ) has to be replaced by \<p n ) and 
the eigenvectors with coefficients c£ by c^. 
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